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Abstract 

Thomas-Fermi theory for Bose condesates in inhomogeneous traps is revisited. 
The phase-space distribution function in the Thomas-Fermi limit is fo(R,p) a 
S(fj, — H c i) where H c i is the classical counterpart of the self-consistent Gross- 
Pitaevskii Hamiltonian. No assumption on the large N-limit is introduced 
and, e.g the kinetic energy is found to be in good agreement with the quantal 
results even for low and intermediate particle numbers N. The attractive case 
yields conclusive results as well. 



1 



I. INTRODUCTION 



The recent discovery of the Bose-Einstein condensation of magnetically trapped atoms 
has spurred a huge amount of theoretical investigations. Most of them are based on the 
Gross-Pitaevskii equation (GPE) |I[ which is the mean-field equation for the condensate 
wave function (order parameter). The experimental conditions are such that the atomic 
gas is at very low density and therefore the mean-field approximation gives indeed excellent 
results Since the number N of atoms involved is generally large, it is natural that also 

the Thomas-Fermi (TF) approximation is applied quite extensively. This has the advantage 
of yielding in most cases explicit analytical results of great physical transparency. However, 
with respect to Fermi statitstics, the TF-approach to Bose-Einstein condensation shows 
some peculiarities which in the past have, in our opinion, not fully been born out. It is the 
purpose of the present paper to further elaborate on the TF-approach to inhomogeneous 
Bose systems. We deliberately restrict ourselves here to the TF-limit of the GPE at zero 
temperature. Finite temperature as well as more elaborated theories like the Bogoliubov 
approach may be the subject of future work. 

In detail the paper is orgasnized as follows. In Section 2 the theoretical aspects of the 
TF approximation to the GPE are presented in detail. In Section 3 the numerical results 
obtained with the TF method are compared with the ones coming from the GPE. The last 
Section is devoted to discussions and an outlook. 



II. GROSS-PITAEVSKII EQUATION. THOMAS-FERMI LIMIT 

As mentioned in the introduction, the basic equation for Bose condensed atoms confined 
by magnetic traps is, in the low density limit, given by the GPE for the wave function of 
the condensate: 

(-^-A + V ex + g\ij\ 2 )ij = nij (1) 

where V ex is the external potential which for simplicity we have considered to be a spherical 
harmonic oscillator (for non-spherical geometry see remarks at the end of the paper). The 
coupling constant is given by g = Antfa/m with m the atomic mass and a the scattering 
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length. The chemical potential p is identical with the lowest eigenvalue of the self-consistent 
potential 

V = V ex + g\^\ 2 . (2) 
It is useful to note that eq.(|J) can also be rewritten as an equation for the density p = \ip\ 2 : 



2- 

P 2 P 



+ V ex + gp - p = (3) 



2m 4 

In the large N limit one can drop in (P the kinetic energy or, equivalently, in (0) the 
gradients terms of the density. This leads to p = (p — V ex )/g what is known in literature 
|3|-^|] as the TF solution of the GPE. However, for moderate particle numbers the kinetic 
energy is not negligeable and there is no reason for dropping it in the TF limit. Also in the 
attractive case the kinetic energy is of crucial importance to avoid collapse. 

As it is well known from the case of the Fermi systems ||, the TF-approximation is based 
on the assumption of a slowly varying potential so that its gradients can be neglected to 
lowest order. The TF approximation of the density matrix corresponding to a single wave 
function is, as a matter of fact, well known Q and given by: 

lim mr)iP(r')} w = NcS(p - H c ) + 0(h 2 ) (4) 

where the index " W" stands for the Wigner transform |J and H c is the classical hamiltonian: 
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H c = + V ex (R) + gp{R)- The constant c{p) of dimension of an inverse of energy is 
determined from the condition that the wave function must be normalized leading to : 

1 r dRdp 



c J {Inn) 6 

As a consequence the chemical potential p must be determined from an independent quanti- 
zation condition (see below). We will, however, see that in the limit iV — > oo the expression 
for p coincides with the usual. 

It is clear from (|^) that the kinetic energy has been properly included (as usual in TF 
theory). One way to understand expression (f|) is to write eq.(P in the form (p — H)p = 
with p = \ip){ip | the density matrix. Wigner transforming this equation and remembering 
that to lowest order in % the Wigner transform of a product of operators is the product of 
their Wigner transforms ||, one arrives with x8(x) = at eq.(§). 



A. Self-consistent solution. Repulsive case (g > 0) 



Let us first consider the solution of the self-consitent problem at the TF level defined by 
eq.(^) for the repulsive case i.e g > 0. From (Q) we obtain for the density: 

. „, f dp „ , „ . Ncm , 

" (i{) = /(W /o(i{ ' p) = ^ Po(ii) (6) 

where the local momentum is given by : 



Po (R) = pmin-Vex-gp) (7) 
The self-consistency between eqs.(||,0) is easy to solve analytically and we obtain : 



where 

It is to be noted that p(R) is defined only within the classical region limited by p — V ex (R) = 
0. It is straightforward to expand p in the repulsive case g > for large values of N : 

P"lb-V„)-j^(M-V„)> + (10) 

It is satisfying to see that to leading order one recovers the result corresponding to the total 
neglect of kinetic energy in (P (see introduction). 

The normalization is directly determined from eq. @ 

- = ^3 / dRpoR) (11) 
c 2n 2 h J 

or equivalent ly 

N = J dRp(R) (12) 
Explicitly one obtains from this equation: 



mu 2 J I 48 r 64 



Using in (|i~2[ ) the asymptotic expansion (|T^) yields 



^mu?) TE~gN\ ~ 7 Kg 2 N 2 ) (H) 



which also can directly be derived in expanding (|13f) . 

From Eqs. (p~3| ) or (|14D we can determine the normalization constant c as a function of \x. 
From (14) we see that to lowest order the normalization constant c drops out and thus in 



this limit the chemical potential is, as usual, determined by the particle number condition. 
However, as we will see, via the quantization condition [i depends on c and thus we can 
consider (|T3"| ) as determining the normalization in any case. 

B. Chemical potential and quantization 

The semiclassical density matrix (|j) corresponds to a single wave function. In such a 
case the energy must be determined independently from a quantization condition. This is in 
fact well known ||. Formally in our case this is necessary, since we have two open constants 
fi and K. The equation for the chemical potential \x = dE/dN is also of no help, since it 
is equivalent to the particle number condition (jTJD The standard semiclassical quantization 
procedure is given by the WKB method. However, in order to have a more explicit formula, 
we here also apply a slightly simpler method, aplicable to the lowest state in a single particle 
potential 0. To this end we calculate the smooth accumulated level density (number of 
states) in TF approximation: 

= J (15) 

For a spherical harmonic oscillator (H.O.) this gives: 
Taking for E in (|16|) the H.O. eigenvalues 



E -> E K = (K+ ^)hu (17) 

with K = 2n + I and inserting (|1(],|T7]) in the left hand side of (|T5| ) yields a semiclassical 
quantization rule which becomes exact in the 3-D spherical H.O. case. It represents an 



approximate quantization relation for an arbitrary potential where the quantized energies 
very well reproduce the centroid of major shells. This has been tested numerically on a 
potential of Woods-Saxon type of nuclear dimensions [0] . It is evident that in the 1-D case 
the same procedure leads to the exact WKB quantization rule. For the 3-D case this modified 
quantization prescription is slightly less accurate than WKB for the lowest eigenvalue but 
has the advantage to be easier and to be readly applicable also to the deformed case ||. In 
the present problem the eigenvalue \i is then determined by: 
27 1 



^3 J dR Po 3 (R) 

32 f 27 f Kg 2 N 2 \^ 2 5/2 19 ( Kg*N* \W 3/2 7 / KgW \W 1/2 



7T 

12 f 27//\/r'.V-\ ; : 

irW\~ 80 V 4 ) fi ' ~24{ 4 ) ^ ' ' 16 V 4 1 '"' 
, 3 Kg 2 N 2 ( Kg 2 N 2 \ 2 . / g 

+ s — 1 — I — 1 hA "J arcsm Vl2^ 

,2 JU2 x 3 



1 ( Kg 2 N 2 \" . I fi \ 

H h A* arcsm . / 2 .. 2 > 

16 V 4 'V \ *s^! + w J 



4 + M' 

where we have used Po{R) from (0) with (H). To second order we obtain from (|i~8l): 



27 4 / 2m 1 \ 3 / 2 \ 3 / 2 9/2 16 / 48 fi \ 
8 vrft 3 VV K gNj Kmou 2 ) ^ 315 1 11 Kg 2 N 2 ) 1 j 

To leading order in the large N limit we obtain from ( |T3"| , |T4] ) and ( JT5|JT~D| ) : 

15\ 2 / 5 /mw 2 \ 3 / 5 / _\ 2/5 



{—) r 

K <> - Graf w"° 3 < 20 > 



This completes the solution in the large N limit. 



For later comparison let us also give the standard WKB quantization rule [10] which we 
want to evaluate to leading order : 

h 2 1/4] 



7i 1 2m r 2 , 
, , dr 

2 V ft 



H - V{r) 



2m r 2 



(21) 



with (see eq.fliHD) 



V~V(r) = 1 ^- 2 ( f i-V ex (r)) 2 (22) 
The classical turning points r\ and r2 are determined from the solution of the cubic equation: 



y/Kg*N 2 K ex ' 2V2mr 
With d|) and (H) eq.© can be solved for fi. 



1 1 ft 2 1 

(M-^0 = -o\ 77-- (23) 
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C. Kinetic energy 



One of the main difficulties with the standard iV — ■> oo TF limit treated in the literature 
(see introduction) consists in the inability to calculate the kinetic energy ||. In our approach 
this does not cause any particular problem and one directly obtains 



E 



kin 



For example to lowest order we obtain from 



PO = rg= — (Ho ~ V ex )Q(fl - Vex) (25) 

VK gN 



what yields for the kinetic energy per particle, using (£(J 



41 = 272^ 27 (2g) 

N 32 m V 2m 32 1 ; 

This simple result must be contrasted with the usual statment that in the N — > oo TF limit 
the kinetic energy can not be evaluated, since it diverges ||. 
Indeed one can write the kinetic energy as: 

E km = ^-JdR\Vij\ 2 (27) 

In the N oo limit we have: 



IpN^oo = J -(fid - 14b) (2£ 

V ^ 



and one can readily verify that with fl28|) of of (|27D diverges logaritmically. This result 
obviously is in contradiction with ( |26D and we shortly want to elucidate the underlying 
reason. To this end we first rewrite (E7I) in a different but obviously equivalent way: 

Eun ^j dRm ^J^^MR, P ) (29) 



with /o given by the Wigner transform of the density matrix corresponding to (|28 



fo(R, p)= J dse-'P^N^R + -^^(R - |) (30) 

Since Jo ^ fa we argue that is not the correct ft — > limit of the distribution function 
because it is not solution of the h — > limit of the Schrodinger equation ([]]) 



(H c -ti)f o = (31) 
Only (^) is the correct solution of this equation which yields for large N: 

A ~ = '" , (ji^ ta - ,u, -s) (32) 

One checks that upon projection onto r-space fl3~2|) gives the correct lowest order expression 
for the density (see eq.( fL0|) ). Therefore both Wigner functions (j30| ) and (|32|) yield the same 
leading order density. However, in spite of being a very suggestive non-local generalization 
of the lowest order local density expression, eq.(|3~0"D has to be rejected on the above given 
grounds and the divergency of (p 7 ^ ) is an artifact. On the contrary the lowest order contri- 
bution to the kinetic energy is given by (PB|). Via (|T^,|T^) it is straightforward to calculate 
the next to leading order correction to (|26|) . It should, however, be remembered that 1/N 
correction do not go in parallel with powers in H and that 1/N corrections also can come 
from h 2 corrections to (f|) which involve second order gradients of the potential. In any case 
the Wigner-Kirkwood expansion of the density matrix is an asymptotic expansion which in 
no way can recover the nonanalytic behavior in h of the quantal solution. In the present 
problem the nonanalyticity in H of the quantal solution entails a nonanalytic behavior in 1/N 
(see eq.(16) of ref. 0) and therefore a WK-expansion can never recover the quantal behavior 
in 1/N. It is well known that an asymptotic expansion has to be stopped at a point where 
the difference to the exact solution is minimal. Afterwards the expansion starts to diverge 
again. In this work we do not intend to develope a systematic expansion simultaneously in 
h and 1/N. We rather want to give a complete solution to lowest order in H, i.e. on the TF 
level. 



D. The attractive case (g < 0) 

It seems that recently, Bose-Einstein condensation has been observed also for the case of 
negative scattering length ( u Li atoms ) JTIJ]. 

For g < the Gross-Pitaevskii approach leads to metastability for particle numbers 
N < 1400 0. For large particle numbers the system collapses. For the attractive case 
(a < 0) the correct treatment of the kinetic energy is crucial in the TF-limit, since otherwise 
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no stability can be achieved. Formally the TF solution for the density is the same as in @ 
with, however, the sign of the first member reversed : 



Contrary to the repulsive case no large N expansion is possible here. Therefore the TF 
solution has to be considered in full. In the next section comparison with quantal results 
will be given. 



In this section we proceed to a detailed numerical comparison of the semiclassical ap- 
proximations with the exact quantum mechanical results. Along this section energies and 



lengths are given in harmonic oscillator units : Tiu and duo — \Jh/2mu, respectively. First 
in Table 1 we present the chemical potential (/x) and the kinetic (ejy„), harmonic oscillator 
(cho) an d the self-interaction (e p ) energies per particle calculated quantally and in the full 
Thomas- Fermi approximation (p|,|l3l|l8|) as a function of the number of atoms enclosed in the 
trap. We have considered Cs atoms (as was done in Ref. [0), the frequency of the harmonic 
oscillator has been chosen to be u = 207rs~ 1 and the scattering length to be a = 3.2xl0 _9 m. 

In Table 2 we present the results for the chemical potential and the kinetic energy per 
particle number beyond 20000. In addition to the quantum mechanical and the full Thomas- 
Fermi results, we also include the results for the large iV limit ( PB|) and those obtained using 
the WKB quantization rule (in the large N limit). Notice that in the large N limit, the 
WKB chemical potential coincides with the TF one and the kinetic energy is given also by 
( p6|) but with c replaced by the one calculated from K via eq.fl2"T|). From Tables 1, 2 we see 
that for instance the results of the full TF solution are in very satisfaying agreement with 
the quantal results over the whole range of particle numbers. 

The numbers presented in Table 2 indicate that the asymptotic values of the chemical 
potential and the kinetic energy are obtained only for very large number of particles (N ~ 
10 5 — 10 6 ). We also realise that the WKB quantization rule yields quite similar results though 
in fact slightly worse ones than our simpler TF-quantization rule ( JT8|JT9| ) . Though globally 
the semiclassical results of Tables 1 and 2 are quite satisfactory, one nevertheless remarks 




(33) 



III. NUMERICAL RESULTS 
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some unexpected features. For instance the kinetic energy in the TF approximation is larger 
than the exact values for small numbers of particles whereas it undershoots the quantum 
values quite considerably in the large N limit. We will come back to a more detailed analysis 
of this behavior in the discussion section. 

Next let us compare in Fig.l and 2 the densities (normalized to unity) in TF- 
approximation and calculated exactly for small (200) and large (200000) particle numbers. 
As expected, the TF densities almost agree with the quantal ones for very large particle num- 
bers. In view of the still quite reasonable expectation values shown in Table 1 for N=200, 
the strong deviation of the TF density from the quantal result is somewhat a surprise. How- 
ever, one always should remember that the TF-solution for the densities is to be understood 
as a distribution (see eq|4|) which for expectation values of "slowly varying" operators can 
still yield very reasonable values in spite of the fact that the detailed shape may only be a 
charicature of the exact one. 

In Figs. 3 and 4 we show the self-consistent potentials V = V ex + gp corresponding to the 
densities of Figs.l and 2. Not astonishingly V deviates from the harmonic oscillator only 
slightly for N = 200, both quantally and semiclassically. On the contrary for iV = 200000, 
the potential V deviates strongly from V ex being practically a constant equal to /x up to 
the classical turning point from where the harmonic oscillator takes over quite abruptly. 
Again, both quantal and TF-solution are in close agreement. Fig.4 also teaches us why 
the TF-approximation @ to the quantal distribution function is very good for large N. The 
distribution function corresponds to a wavefunction with very large energy \x. In phase space 
it therefore is very much concentrated around the surface of the hypersphere with radius \l. 

Let us now present the attractive case for the same atoms and external potential with 
, however, the scattering length a = — 1.0xl0 _9 m. In Table 3 we again show chemical poten- 
tial and kinetic, harmonic oscillator and self-interaction energies per particle as a function of 
the particle number in TF and quantal calculation. For small particle numbers (N < 1000) 
the agreement of TF with the quantal case is of similar quality as in the repulsive case. How- 
ever for N > 1000 the agreement quickly deteriorates, indicating that the whole mean field 
approximation breaks down. Indeed even quantally the solution of the GPE ([!]) becomes 
unstable for N > 1500 for Cs atoms. 
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In Figs. 5, 6, 7 and 8 we also show the densities (normalized to unity) and self-consistent 
potentials for the particle numbers iV = 250 and iV = 1500. We see that, whereas the case 
iV = 250 is not dissimilar to the corresponding one with a > 0, for N = 1500 the situation 
becomes quite unfavorable for the TF approximation. This is for instance manifest in looking 
at the graph for the densities. In the attractive case TF and quantal solutions diverge with 
increasing N whereas in the repulsive case they converge. 

In Figs. 9 and 10 we plot the kinetic energy density per particle (t/N) calculated quantally 
and in the TF-approximation for iV = 200000 in the repulsive case and for iV = 250 in the 
attractive case. In these Figures the quantal kinetic energy density is given by: 



r=\^\ 2 - 1 -Ap ==I 



(Vp) 



2 



Ap 



(34) 



P 

in order to compare with the TF one according to Ref. 

For large number of particles, when the density profile has a relatively flat region at the 
interior (see Fig. 2), the quantal kinetic energy density is peaked at the surface (see eq.(|34"|)) 
whereas the TF one is rather a bulk term (see eq. (|24|) ) . Inspite of the rather different form 
of the quantal and TF kinetic energy densities in this case, the corresponding integrals are 
in good agreement (see Tables II). This fact points again to the distribution character of 
the TF-kinetic energy density. However, if the number of particles is small, the quantal and 
TF kinetic energy density profiles are quite similar. This is due to the fact that the particle 
density fall-off abruptly from R = (see Fig. 5) and consequently its derivatives contribute 
in all the range of R. 



IV. DISCUSSION AND OUTLOOK 



In the preceding sections we have derived the Thomas-Fermi approximation i.e. the 
h — > limit of the density matrix corresponding to the wavefunction of the Bose condensate 
of atoms confined by magnetic traps. We have pointed out some misconceptions on this 
point which appeared in the past in the literature which for instance prevented the direct 
calculation of the kinetic energy in the large N limit. On the contrary with our Thomas- 
Fermi approach the evaluation of the kinetic energy causes no problem and the results are 
globally in quite satisfactory agreement with the quantal solution of the Gross-Pitaevskii 
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equation. However from Table 2 we see that the kinetic energy in the TF limit does not 
have the correct asymptotic behavior as a function of particle number. As was shown by 
Pitaevskii and Stringari |5), this is due to nonanalytic (logarithmic) corrections which can 
not be accounted for by a pure TF approach and needs a partial resummation of all orders 
in H. However, for particle numbers where the kinetic energy represents a significant fraction 
of the total energy the TF expression yields very satisfying results for E^ n . This example 
shows again that the semiclassical approximations are a powerful tool but not devoid of 
subtleties and pitfalls. As a matter of fact also in this paper we, for simplicity, avoided 
to develope the full complexity of the theory. One major simplification resides in the fact 
that we assume a spherical trap. This results in an isotropic momentum distribution fo a 

2 

S(/i — H c ) where H c = ^— + V is the classical Hamiltonian. Deforming the trap leads to a non 
trivial modification of the theory, since squeezing the condensate wavefuntion in one direction 
and relaxing in the other entails in turn a deformation of the momentum distribution which 
is opposite to the spatial one, i.e. momenta are strongest in the squeezed direction and 
lowest in the long direction of the deformation Our TF approach can also be useful 

for the evaluation of collective excitations of droplets of small or intermediate sizes. Such 
a situation is in fact well known from the zero sound giant quadrupole vibrations of finite 
nuclei where to first approximation the restoring force of the vibration is given by the energy 
stored in the deformation of the Fermi sphere || . On the other hand the deformation of the 
momentum distribution of condensed atoms has also been revealed experimentally in directly 
measuring the momentum distribution of the expanding particles, once the deformed trap 
has been turned off. The detailed determination of the anisotropy of the momenta (which 
may be position dependent) is theoretically a not completely trivial task in the general case 
and we will elaborate on this in future work. In the present case, however, there exists 
an evident first guess of the momentum deformation which results from a scaling argument 
of the harmonic oscillator coordinates. Assuming a prolate quadrupole deformation in the 
^-direction we have to replace the classical hamiltonian in (^) by : 

2 2 n 

^0 / 2 , 2\ , ^0 2 



H c = 

2m 



( Px 2 +Py 2 ) + ^pA + V(^x, ^y, ^z) (35) 



where the ratios — and — are the frequency relations in z and x, y with respect to the 

LJq (jJq 

spherical case (wo)- From (^) one easily calculates the so-called aspect ratio in the TF- 
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approximation 



% = ^ (36) 

Px UJ Z 



a result which has been given previously 0. One other important consequence of the mo- 
mentum deformation is that with (|35|) the moment of inertia of the condensate becomes equal 



to the irrotational flow value ||14|| . On the contrary using ([|) with the isotropic momentum 
distribution the rigid momentum of inertia results. Consequently the deformed case needs 
more detailed studies which we reserve to future work. It is also evident that the present 
TF approach can be extended to finite temperature and to the Bogoliubov theory. 

Another interesting subject of a more formal aspect is the evaluation of the ^-correction 
to the present lowest order theory. In principle this can easily be performed in posing in (|3|) 
P = Po + ^ 2 P2 and = /io + ^ 2 /i2 and properly sorting out different powers in h. However 
the proper elimination of divergencies and handling the normalization ( fl2|) and quantization 
(|18|) are slightly subtle problems. 

Investigations on the above mentioned directions are in progress. 
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TABLES 



TABLE I. Chemical potential (/i), kinetic energy (e^A, harmonic oscillator energy (eno) 
and self-interaction energy (e p ) per particle in harmonic oscillator units (hu) calculated 
quantally (QM) and in the Thomas-Fermi approaximation (TF) for several numbers of atoms 
in the traps. The frequency of the harmonic oscillator is uo= 20 n s _1 and the scattering 
length is a= 3.2xl0~ 9 m. 

TABLE II. Chemical potential (fx) and kinetic energy (e kin ) per particle for large number 
of atoms in the trap (N). The chemical potential is calculated quantally (QM), with the full 
Thomas- Fermi approach (TF) and with the asymptotic formula for large N (TFjv^oo)- The 
kinetic energy is obtained quantally, with the exact TF approximation, with the asymptotic 
TF for large N and using the WKB quantization in the limit of large number of atoms. The 
frequency of the harmonic oscillator and the scattering lenght are the same as in Table I. 

TABLE III. The same as in Table I but with a scattering length a= -l.OxlCT 9 m. 
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TABLE I 



N m(QM) ^(TF) e kin (QM) e kin (TF) e HO (QM) e HO (TF) e p (QM) e p (TF) 



200 


1.688 


1.642 


0.696 


0.700 





.811 





.804 


0.080 


0.069 


400 


1.806 


1.766 


0.654 


0.661 





.865 





.851 


0.144 


0.127 


600 


1.927 


1.877 


0.622 


0.630 





.912 





.894 


0.196 


0.177 


800 


2.036 


1.978 


0.597 


0.603 





.955 





.934 


0.242 


0.220 


1000 


2.134 


2.071 


0.575 


0.581 


0. 


.944 





.970 


0.282 


0.260 


1200 


2.225 


2.157 


0.557 


0.561 


1 


,031 


1. 


.005 


0.319 


0.296 


1400 


2.310 


2.238 


0.541 


0.544 


1 


,065 


1 


,037 


0.352 


0.329 


1600 


2.389 


2.315 


0.528 


0.528 


1 


.065 


1 


.068 


0.382 


0.359 


1800 


2.464 


2.388 


0.515 


0.514 


1 


.127 


1 


,097 


0.411 


0.388 


2000 


2.535 


2.457 


0.503 


0.502 


1 


.158 


1 


.125 


0.437 


0.415 


4000 


3.112 


3.025 


0.431 


0.417 


1 


.395 


1 


.356 


0.643 


0.626 


6000 


3.550 


3.461 


0.390 


0.369 


1 


.577 


1 


.536 


0.792 


0.778 


8000 


3.914 


3.825 


0.363 


0.336 


1 


,729 


1. 


.687 


0.911 


0.901 


10000 


4.231 


4.142 


0.343 


0.312 


1 


,862 


1 


,820 


1.013 


1.005 


12000 


4.513 


4.426 


0.327 


0.293 


1. 


,981 


1 


.939 


1.103 


1.097 


14000 


4.770 


4.684 


0.314 


0.277 


2 


,089 


2 


,047 


1.184 


1.180 


16000 


5.007 


4.921 


0.303 


0.265 


2 


.189 


2, 


.147 


1.258 


1.255 


18000 


5.228 


5.143 


0.294 


0.254 


2 


.282 


2 


.240 


1.326 


1.324 


20000 


5.435 


5.350 


0.285 


0.244 


2 


.369 


2 


.328 


1.390 


1.389 
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TABLE II 



N fi(QM) //(TF) ^(TF^oo) e kin (QM) e kin (TF) 

e fcin(TF7v^oo) e fein(WKB7v^oo) 



20000 


5.435 5.350 


5.196 


0.285 


0.244 




0.256 




0.238 


30000 


6.322 6.242 


6.111 


0.255 


0.210 




0.218 




0.202 


40000 


7.051 6.973 


6.856 


0.236 


0.187 




0.194 




0.180 


50000 


7.677 7.603 


7.496 


0.222 


0.173 




0.177 




0.165 


100000 


10.231 9.972 


9.891 


0.182 


0.133 




0.134 




0.125 


150000 


11.763 11.701 


11.633 


0.162 


0.113 




0.114 




0.106 


200000 


13.170 13.112 


13.051 


0.149 


0.101 




0.102 




0.095 


250000 


14.381 14.326 


14.270 


0.140 


0.093 




0.093 




0.087 








TABLE III 












QM 


TF 


QM TF 


QM 


TF 


QM 


TF 






TV 250 


250 


500 500 


1000 


1000 


1500 


1500 






H 1.424 


1.437 


1.338 1.369 


1.120 


1.212 


0.691 


1.009 






e k in 0.788 


0.774 


0.815 0.802 


0.926 


0.872 


1.240 


0.974 






e HO 0.723 


0.726 


0.691 0.702 


0.613 


0.646 


0.472 


0.579 






e p -0.039 


-0.032 


-0.084 -0.070 


-0.209 


-0.153 


-0.511 


-0.272 
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FIGURE CAPTIONS 



Figure 1. Density (nomalized to unity) of 200 atoms in a spherical trap in aj^ units) as 
a function of the distance (in a^o units) in the repulsive case calculated from the solution 
of the GPE (solid line) and using the TF approach described in the text (dashed line). 

Figure 2. The same as Figure 1 but with 200000 atoms in the trap. 

Figure 3. Self-consisten potential (in ftuj units) corresponding to a spherical trap con- 
taining 200 atoms as a function of the distance (ano units) in the repulsive case calculated 
from the solution of the GPE (solid line) and using the TF approach described in the text 
(dashed line). 

Figure 4. The same as Figure 3 but with 200000 atoms in the trap. 

Figure 5. The same as Figure 1 but with 250 atoms in the trap in the attractive case. 

Figure 6. The same as Figure 1 but with 1500 atoms in the trap in the attractive case. 

Figure 7. The same as Figure 3 but with 250 atoms in the trap in the attractive case. 

Figure 8. The same as Figure 3 but with 1500 atoms in the trap in the attractive case. 

Figure 9. Kinetic energy density per particle of a spherical trap (in hu a]j 3 units) 
containing 200000 atoms as a function of the distance (in clho units) in the repulsive case 
calculated from the solution of the GPE (solid line) and using the TF approach described 
in the text (dashed line). 

Figure 10. The same as Figure 9 but with 250 atoms in the trap in the attractive case. 
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